Salicylic acid carboxyl glucosyltransferase UGT87E7 regulates disease resistance in Camellia sinensis

Abstract Plant immune response following pathogenic infection is regulated by plant hormones, and salicylic acid (SA) and its sugar conjugates play important roles in establishing basal resistance. Here, the important pathogen Pseudopestalotiopsis camelliae-sinensis (Pcs) was isolated from tea gray blight, one of the most destructive diseases in tea plantations. Transcriptomic analysis led to the discovery of the putative Camellia sinensis UDP-glucosyltransferase CsUGT87E7 whose expression was significantly induced by SA application and Pcs infection. Recombinant CsUGT87E7 glucosylates SA with a Km value of 12 µM to form SA glucose ester (SGE). Downregulation reduced the accumulation of SGE, and CsUGT87E7-silenced tea plants exhibited greater susceptibility to pathogen infection than control plants. Similarly, CsUGT87E7-silenced tea leaves accumulated significantly less SA after infection and showed reduced expression of pathogenesis-related genes. These results suggest that CsUGT87E7 is an SA carboxyl glucosyltransferase that plays a positive role in plant disease resistance by modulating SA homeostasis through a mechanism distinct from that described in Arabidopsis (Arabidopsis thaliana). This study provides insight into the mechanisms of SA metabolism and highlights the role of SGE in the modulation of plant disease resistance.


Introduction
Plants face various environmental pressures while growing in nature, including biotic and abiotic stresses. Over the long course of plant adaptation to the environment, plants have acquired a multistage immune system to protect themselves against pathogens. Plant immune response includes pathogen-associated molecular pattern-triggered immunity, which allows plants to fend off a large number of potential pathogens (Li et al., 2005), and effector-triggered immunity (Jones and Dangl, 2006), which relies on proteins encoded by resistance (R) genes to induce a stronger and more durable immune response (Chisholm et al., 2006). When plants are exposed to pathogens, some responses are limited to the infested damaged organ, while other responses systemically spread far from the infested organ and affect the whole plant (Conrath et al., 2015;Vlot et al., 2021). These latter responses include systemic-acquired resistance (SAR) and induced systemic resistance (ISR), depending on the site of induction and the lifestyle of the pathogenic microorganism (Choudhary et al., 2007;Vlot et al., 2021). Small metabolites, such as salicylic acid (SA; Lim et al.,2020), glycerol-3-phosphate (Chanda et al., 2011), azelaic acid (Jung et al., 2009), dehydroabietinal (Chaturvedi et al., 2008), pipecolic acid (Pip; Návarová et al., 2012) and N-hydroxy-pipecolic acid (NHP; Chen et al., 2018;Hartmann et al., 2018) have been reported to be involved in long-distance communication and plant immunity (Shah et al., 2014;Vlot et al., 2021). Among them, SA plays a vital role in the activation of defense mechanisms during pathogen infection (Vlot et al., 2009) and is associated with the accumulation of pathogenesis-related (PR) proteins (Ross, 1961;Durrant and Dong, 2004).
The plant hormone SA plays a regulatory role in many physiological and biochemical processes, serves as a critical signal for the activation of disease resistance in plants, and participates in plant immune response (Vlot et al., 2009;Ding and Ding, 2020). Once synthesized, free SA may undergo several biologically relevant chemical modifications, including glucosylation and methylation. Most modifications render SA inactive, permitting the fine-tuning of its accumulation, function, and/or mobility (Dempsey et al., 2011). In plants, small hydrophilic molecules such as glucose are often conjugated to SA to aid in transport and storage because of their reactivity and hydrophilicity (Thompson et al., 2017). MeSA, the methylated derivative of SA, has been reported to act as a signaling compound in SAR in Nicotiana benthamiana (Park et al., 2007). A different study, however, indicated that neither MeSA nor JA is essential for systemic immunity, emphasizing the crucial role of SA in Arabidopsis (Arabidopsis thaliana) (Attaran et al., 2009). SA glucose exists in two forms, SA 2-O-b-D-glucoside (SAG) and SA glucose ester (SGE). In Arabidopsis, in vitro catalytic analyses determined that the glycosyltransferase UGT74F1 can form SAG, whereas UGT74F2 produces both SAG and SGE (Lim et al., 2002). The A. thaliana glucosyltransferase UGT76B1mediated glucosylation controls the levels of active SA, NHP, and isoleucic acid in concert to balance plant defense and growth Cai et al., 2021;Hõrak, 2021;Mohnike et al., 2021, Zeier, 2021. SAG is thought to be a more stable storage form for small phenolics, and its levels have been shown to increase in parallel with free SA in plant defense responses, suggesting that it is associated with plant defense (Enyedi et al., 1992;Chen et al., 1995). SGE is thought to be a high-energy compound and biosynthetic intermediate (Bowles et al., 2006), but the formation of SGE and its role(s) in plant defense is still not understood (Boachon et al., 2014).
Tea is considered to be the most popular nonalcoholic beverage and is produced from leaves of the tea plant (Camellia sinensis) . With the continuous expansion of tea plantation area globally, damage from tea diseases has seriously affected tea production and quality. The leaf disease, tea gray blight, also known as tea shoot wilt, is one of the main diseases of tea plants in China . Its common pathogen has also reduced the production of apple and guava, causing considerable economic losses (Harteveld et al., 2014;Solarte et al., 2018). Tea gray blight is instigated by Pestalotiopsis-like species; the infection appears initially as small brown marks on the leaves and develops into large round brown necrotic lesions, severely affecting the quality of tea production.
In this study, Pseudopestalotiopsis camelliae-sinensis (Pcs) was isolated from tea gray blight, and its cytological characteristics were investigated before pathogenicity tests. Transcriptomic analysis was performed on Pcs-infected tea plants, the C. sinensis UDP-glucosyltransferase (UGT) CsUGT87E7 was identified, and its expression was shown to be strongly induced by infection. In vitro assays showed that CsUGT87E7 catalyzes the glucosylation of SA, and its main product was identified as SGE. Downregulation of CsUGT87E7 expression reduced the accumulation of SGE and increased the susceptibility of tea plants to pathogen infection. After infection, CsUGT87E7-silenced tea leaves accumulated significantly less SA and showed lower expression of PR genes. These results suggest that CsUGT87E7 plays a key role in plant disease resistance through a mechanism distinct from that previously described in Arabidopsis. This study provides insight into SA metabolism and highlights the role of CsUGT87E7 in the formation of SGE and the modulation of plant defense.

Pathogen isolation and pathogenicity analysis
Tea plants with typical gray blight disease symptoms (concentric gray rings with scattered black granules) under field conditions were selected for pathogen isolation and identification (Figure 1, A). The pathogen was purified from the diseased tea leaves and identified to the species level based on its morphological (Figure 1, B) and molecular characteristics. Genomic DNA was extracted from the purified pathogen and amplified, and the internal transcribed spacer (ITS) gene region was sequenced using primers ITS1 and ITS4 (White et al., 1990). The ITS sequence of the pathogen showed 99% similarity with Pcs sequences. We therefore confirmed that the fungus was Pcs based on its morphological and molecular characteristics.
The pathogenicity of the Pcs isolate was confirmed using Koch's postulates. The leaves of healthy tea plants developed visible gray blight symptoms after inoculation with the Pcs strain. Five to ten days after inoculation, disease symptoms were noted on the epidermal surface, including extensive cell necrosis and lesion development (Figure 1, C). These symptoms were consistent with the symptoms observed in the above-mentioned gray blight disease under field conditions (Figure 1, A). Therefore, Koch's postulates confirmed that the pathogen re-isolated from the diseased tea leaves was the Pcs pathogenic fungus.

Fungal infection triggers the formation of SGE
It has been proposed that SA acts as an endogenous signal responsible for inducing SAR in plants (Gaffney et al., 1993). We therefore measured SA content by liquid chromatography-mass spectrometry (LC-MS) in tea plants after Pcs pathogen infection. The SA content increased continuously in infected plants (about 1.7-to 2.6-fold) compared with control plants that were inoculated with water ( Figure 1, D).
Endogenous or exogenous SA is metabolized to a glucose conjugate in plants ( Lee and Raskin, 1998). To investigate the formation and metabolism of SA after Pcs infection, we measured levels of the SA glucose conjugates SAG and SGE. The SAG content of infected tea plants remained at A B E C F D Figure 1 Isolation of pathogenic fungus from diseased leaves and the levels of SA and two glucosylated SA forms in tea plants after fungal infection. A, Gray blight disease of tea plants in the field. B, The plant pathogenic fungus Pcs isolated from diseased tea leaves and grown on potato dextrose agar (PDA). Photos on the right show conidia and fungal hyphae. Bar = 20 mm. C, Pathogenicity test of fungi on tea plant leaves. D, The relative content of SA after Pcs infection. E and F, The relative contents of SAG and SGESGE after Pcs infection. Control, inoculated with pure water. Data are presented as mean ± SD of at least three replicates. Significant differences between the treatment and the control group were calculated by one-way ANOVA in SPSS 21.0. *P 5 0.05, **P 5 0.01, and ***P 5 0.001. relatively stable levels within 24 hpi and its content significantly increased at 48 hpi compared with that in the control plants (Figure 1, E). SGE was almost undetectable in the healthy tea plants, but it was significantly induced at all selected time points by pathogenic fungus infection (Figure 1, F). Taken together, our results showed that SA is mainly metabolized and stored as SAG in healthy tea plants, whereas SGE formation was influenced earlier than SAG in response to fungal infection. These results suggest that SGE might play some role(s) in the defense against pathogen invasion and motivated us to investigate its formation and function in response to Pcs infection in tea plants.

CsUGT87E7 catalyzes the formation of SGE in vitro
The transcription of glycosyltransferase genes is typically induced by their substrates (Huang et al., 2018;Jing et al., 2019;Zhao et al., 2020). To search for UGTs associated with SGE formation, we performed transcriptome sequencing from 0 to 48 h after the application of exogenous SA. Tea plants sprayed with water served as controls. Transcript levels of putative UGTs were analyzed and 13 UGTs whose expression was substantially induced by SA application were selected as candidates ( Figure 2, A). The full-length sequences of these candidate UGTs were obtained from young leaves of C. sinensis cv. Shuchazao. All UGTs were successfully expressed in Escherichia coli BL21, affinity purified, and used for subsequent enzymatic activity analysis (Supplemental Table S1). The products were analyzed by LC-MS.
The UGT encoded by Tea006702 catalyzed the formation of SA glycoside, as identified by LC-MS. In negative ionization mode, the dominant ion peaks at m/z 137.02 (M-H + -Glc) and 299.01 (M-H + ) were compared with authentic SGE and SAG standards ( Figure 2, B and C). The results showed that the protein encoded by Tea006702 catalyzed the formation of SGE, and it was assigned the name CsUGT87E7 based on the conventions of the UGT Nomenclature Committee (Mackenzie et al., 1997).

Substrate specificity and kinetic constants of CsUGT87E7
The enzymatic activity of CsUGT87E7 was further analyzed by a UDP-GLO glycosyltransferase assay using a broad range of chemical classes, including phytohormones (ABA, IBA, JA, and IAA), flavonols (kaempferol, quercetin, and naringenin), and carboxylic acids (sorbic acid and benzoic acid). UDPglucose served as the donor substrate. CsUGT87E7 did not catalyze the glucosylation of any of these substrates. Furthermore, sugar donor analysis clearly revealed that CsUGT87E7 preferred UDP-Glc (100%) to UDP-Gal, UDP-GA, and UDP-Rha as sugar donor (Supplemental Figure S2).
To determine its kinetic constants, CsUGT87E assay conditions were established using SA and UDP-glucose as the acceptor and donor substrate, respectively. Kinetic properties were determined for SA in the linear range of the enzymatic reaction (2 mg of protein, 30 min reaction time) at optimized conditions (Supplemental Figures S3-S5), and the apparent K M value of CsUGT87E7 for SA was determined to be 12.22 ± 3.22 mM (Figure 2, D). The low K m value for SA strongly suggests that SA is the in vivo substrate of CsUGT87E7. These results indicate that CsUGT87E7 is one of the enzymes involved in the formation of SGE.

Expression of CsUGT87E7 is induced by exogenous SA and Pcs
To determine whether CsUGT87E7 expression was induced by pathogen attack, tea leaves were exposed to a time course of Pcs and SA treatments, and reverse transcription quantitative PCR (RT-qPCR) was used to analyze the relative CsUGT87E7 expression levels. As shown in Figure 2, E and F, CsUGT87E7 was significantly induced by SA application and Pcs infection. Phylogenetic analysis showed that CsUGT87E7 in tea plants belongs to group J ( Figure 3). These results indicate that CsUGT87E7 is a pathogenresponsive gene and probably plays a role in tea plant defense.

CsUGT87E7 catalyzes the formation of SGE in vivo
To further determine whether CsUGT87E7 was involved in SGE formation in the tea plant, its expression was transiently suppressed in tea leaves by the gene-specific antisense oligodeoxynucleotide suppression (AsODN) strategy, as described in Zhao et al. (2019). CsUGT87E7 expression was significantly lower in tea leaves treated with AsODN_CsUGT87E7 for 12 h than in control leaves treated with sense oligonucleotides (sODN), as measured with two independent primers (Figure 4, A).
Next, the SGE content of tea leaves treated with AsODN_CsUGT87E7 was measured by LC-MS. The accumulation of SGE was markedly lower in CsUGT87E7-silenced tea leaves (P 5 0.01) compared with control leaves (Figure 4, B), suggesting that CsUGT87E7 is involved in the formation of SGE in vivo.
CsUGT87E7 positively modulates disease resistance of tea plants Given that CsUGT87E7 could glucosylate SA to form SGE in vivo, the next question was whether CsUGT87E7 was involved in plant defense. To investigate this possibility, we first suppressed the CsUGT87E7 expression of tea seedlings by injecting 1 mL of 10 mM AsODN-CsUGT87E7 solution; control seedlings were injected with sODN. RT-qPCR analysis showed that the relative expression of CsUGT87E7 could be silenced for at least 4 d (Figure 4, C). The seedlings were then inoculated with Pcs, and their phenotypes were observed for 1 week.
Although there were no differences on the leaf surfaces of tea plants with and without Pcs infection (Supplemental Figure S6), fungal invasion was more severe in the CsUGT87E7-silenced tea plants than in the controls in the first day, as observed by wheat germ agglutinin (WGA) staining (Figure 4, D). Four days after inoculation, disease symptoms became more severe in the CsUGT87E7-silenced plants than in the controls, and The GAPDH was used as the internal reference gene. Data are presented as mean ± SD of at least three biological replicates. Significant differences between the treatment and the control group were calculated by one-way ANOVA in SPSS 21.0. *P 5 0.05, **P 5 0.01. silenced leaves than in control leaves, indicating that CsUGT87E7 positively modulates defense against Pcs in tea plants.

CsUGT87E7 affects SA and reactive oxygen species accumulation in infected tea leaves
In plants, the glucose conjugate of SA can exist as either SAG or SGE ( Figure 5, A). To further verify the role of CsUGT87E7-mediated SGE formation in the tea immune response, we inoculated CsUGT87E7-silenced tea plants (AsODN) and control tea plants with a spore suspension of Pcs and measured SA content ( Figure 5, B). As shown in Figure 5, C, the content of SGE was significantly reduced in CsUGT87E7-silenced plants to around 35%-50% of that in the controls during pathogen infection. The SAG content was not changed and even increased at 24 h post inoculation (hpi) by pathogen ( Figure 5, D), indicating that the metabolite flux of SA was shifted to SAG formation due to the silencing of CsUGT87E7. The reduction of SAG formation after 48 h is probably caused by the reduced SA accumulation ( Figure 5, B). Taken together, it is clearly proven that CsUGT87E7 plays a key role in the formation of SGE.
The SA biosynthesis gene ICS1 (Isochorismate synthase 1), the SA biosynthesis regulatory genes EDS1 (Enhanced disease susceptibility 1) and PAD4 (Phytoalexin-deficient 4), and the PR protein genes PR1 and PR2 (Maleck et al., 2000) were studied in the infected leaves in both CsUGT87E7-silenced and control tea plants after pathogen infection. We first measured the effect of gene suppression at each sampling time to confirm that CsUGT87E7 was effectively silenced in the tea plants ( Figure 6, A). We found that the expression levels of the upstream genes ICSI, EDS1, and PAD4 were higher at 0 (without inoculation) and 12 hpi, whereas their expression was significantly reduced in CsUGT87E7-silenced plants when compared with the controls 24 h after infection ( Figure 6, B), indicating that silencing of CsUGT87E7 affects SA biosynthesis in infected leaves of tea plants during pathogen infection. It should be noted that the expression levels of downstream genes such as PR1 and PR2 were markedly reduced in CsUGT87E7-silenced plants when compared with the controls from the earliest time point after infection ( Figure  6, C). Together, these results suggest that CsUGT87E7-mediated SGE formation influences SA homeostasis and disease resistance during pathogen infection.
Reactive oxygen species (ROS) plays a very important role in plant disease resistance, it accumulates rapidly to kill pathogens (Low and Merida, 1996). SA increases ROS levels by inhibiting the activities of ROS-scavenging enzymes (Chen et al., 1993;Durner and Klessig, 1995). To investigate if CsUGT87E7 affects ROS levels, we directly measured the accumulation of ROS. The generation of superoxide anion (O •-2 ) in CsUGT87E7-silenced tea plants was clearly lower than that in the control plants at all selected time points after inoculation (Figure 7, A). In contrast, the H 2 O 2 content is similar in both CsUGT87E7-silenced and control plants (Figure 7, B), indicating that the expression of CsUGT87E7 might affect the accumulation of O •-2 to activate tea plant immune response.

Discussion
CsUGT87E7 is a SA carboxy glucosyltransferase SA exists in plants in free acid and conjugated forms, metabolism is mainly through glucosylation, and methylation or hydroxylation (Lee et al., 1995;Dempsey et al., 2011). Glycosylation of hydroxylated SA (2,3-DHBA and 2,5-DHBA) and methylated SA (MeSA) has been studied in recent years (Zhang et al., 2017;Huang et al., 2018;Chen et al., 2019). Two enzymes (UGT74F1 and UGT74F2) were identified as SA glucosyltransferases among 90 UGTs in Arabidopsis (Lim et al., 2002). UGT74F1 forms SAG, whereas UGT74F2 forms both SAG and SGE from SA and UDP-glucose (Lim et al., 2002;Dean and Delaney, 2008). Recently, UGT76B1 was identified as a novel SAGT in Arabidopsis and shown to catalyze the conversion of SA to SAG (Noutoshi et al., 2012). Although UGT75B1 forms a glucose ester with the carboxyl group of BA, 3-HBA, 4-HBA, and 3, 4-DHBA, very little activity was observed toward SA (Lim et al., 2002). Until now, a specific SA carboxy glucosyltransferase with high efficiency  Table S2). B, Quantitative analysis of SA glucoside ester after CsUGT87E7 suppression. C, The relative expression level of CsUGT87E7 after CsUGT87E7 suppression after 4 d. D. WGA staining results at 12 and 24 h after Pcs infection. Infected fungi are indicated by green staining. Bar = 50 mm at 12 h; Bar = 500 mm at 24 h. E, Disease symptoms after fungal infection were observed under a stereomicroscope. Bar = 1 mm. F, The average lesion diameter was calculated at 4 dpi. G, The average lesion area was calculated at 7 dpi. Data are presented as mean ± SD of three biological replicates (each with three technical replicates). Significant differences between the treatment and the control group were calculated by one-way ANOVA in SPSS 21.0. **P 5 0.01, *P 5 0.05.
in converting SA to SGE has not been identified, and there have been few studies on SGE function in plants.
Numerous studies have shown that SAGT expression is induced by exogenous SA application or pathogen attack (Kobayashi et al., 2020). Here, we used transcriptomic analysis to identify the putative SA glucosyltransferase CsUGT87E7, whose expression was significantly induced by SA application and Pcs infection (Figure 2, E and F). Both in vitro and in vivo data demonstrate that CsUGT87E7 is a SA carboxy glucosyltransferase (Figures 2, 4). Through phylogenetic analysis (Figure 3), we found that most of the known SA glycosyltransferases, including UGT74F1 and UGT74F2, belong to group L (Lim et al., 2002), although UGT76B1 is found far from them in group H on the phylogenetic tree (Noutoshi et al., 2012). Although phylogenies can be used to predict the functions of unknown genes, some UGTs that cluster together may have very different functions. UGT75B1 and CsUGT75L12 belong to the same group, but UGT75B1 catalyzes the glycosylation of ABA and SA, whereas CsUGT75L12 acts only on flavonoids (Chen et al., 2020;Eudes et al., 2008). CsUGT87E7 identified in the current study belongs to group J (Figure 3), and another member of this group, UGT87A2, has been reported to regulate flowering time and to participate in plant adaptation to abiotic stresses Wang et al., 2012).

SGE formation affects SA homeostasis during pathogen infection
In recent years, SA has been the focus of intensive research because of its function as an endogenous signal and its role in defense responses to pathogens. The accumulation of SA in plants is essential for plant resistance, but most SA in cells is glucosylated and/or methylated (Rivas-San Vicente and Plasencia, 2011). Most endogenous SA, induced by pathogen infection, and exogenous SA application is metabolized into SAG in N. benthamiana (Song, 2006). Like other plant hormones and growth regulators, SAG has been shown to be stored in the vacuole (Dean and Mills, 2004;Dean et al., 2005;Zhang et al., 2017), and it can be converted into SA by a b-glucosidase to serve as a signaling molecule for plant defense response (Seo et al., 1995). In plants, SAG appears to be the predominant metabolite, whereas SGE is a relatively minor metabolite (Dean and Mills, 2004;Song, 2006). Likewise, our results showed that SAG was the predominant glucose conjugate in tea and that SGE was almost undetectable in healthy plants (Figure 1, F). This result is supported by a previous study in which SA was transformed into SAG and methyl salicylate in tea plants (Li et al., 2021).
The existence of a mechanism that releases SA from SAG suggests a possible role for SAG in plant immunity (Hennig et al., 1993). It should be noted that SGE was significantly induced by pathogenic fungal infection, although it was almost undetectable in healthy tea plants (Figure 1, F). This result suggests that SGE plays some role(s) in the response

A B
C D E F Figure 6 The relative expression levels of CsUGT87E7 and PR genes in fungus-infected leaves of control (sODN) and CsUGT87E7-silenced (AsODN) tea plants. A, Relative expression levels of CsUGT87E7 in treated tea plants at different time points. B, Relative expression of SA biosynthesis genes. C, Relative expression of PR1 and PR2. GAPDH was used as the reference gene. Data are presented as the mean ± SD of at least three biological replicates each representing a single infected tea plant. Significant differences between the treatment and the control group were calculated by oneway ANOVA using SPSS 21.0. *P 5 0.05, **P 5 0.01.
to Pcs infection in tea. Silencing CsUGT87E7 in tea plants led to a decrease in SGE. Interestingly, the production of SA was also much lower in CsUGT87E7-silenced plants than in controls, presumably the downregulation of CsUGT87E7 can lead to decreased expression levels of SA-biosynthesis-related genes via feedback inhibition of the structural genes (Yin et al., 2012). In our previous study, downregulation of CsUGT78A14 led to decreased expression of flavonoidrelated genes, ultimately leading to reduced flavonoid content (Zhao et al., 2019). Taken together, our results show that the formation of SGE catalyzed by CsUGT87E7 affects SA homeostasis during pathogen infection, which in turn influences the SA-associated plant resistance in tea plants.
CsUGT87E7 positively regulates disease resistance in tea plants SA-associated plant disease resistance depends largely on the interplay between free and conjugated SA forms (Chivasa and Carr, 1998;Boachon et al., 2014;Ding and Ding, 2020). Generally, plants respond to pathogen infection by rapidly increasing SA content and activating PR genes (Nawrath and Métraux, 1999;Park et al., 2004;Vlot et al., 2021). Uninfected regions of the plant are induced to generate long-term resistance (Vlot et al., 2021). In this study, we showed that CsUGT87E7 was induced by pathogen infection (Figure  2, F). In Arabidopsis, SA glucosyltransferase1 (AtSGT1) was rapidly induced by exogenous SA and infection with a bacterial pathogen, indicating that AtSGT1 expression is an early disease response (Song, 2006). In Arabidopsis, UGT74F1 forms only SAG, whereas UGT74F2 forms both SAG and SGE (Lim et al., 2002). The ugt74f1 mutant showed enhanced disease susceptibility, whereas the ugt74f2 mutant showed enhanced resistance to the same pathogen (Boachon et al., 2014) by negatively influencing the accumulation of free SA. By contrast, the expression levels of PR genes were all reduced after Pcs infection in CsUGT87E7-silenced tea plants, leading to more severe disease symptoms (Figure 4). These results indicate that CsUGT87E7 positively regulates disease resistance in infected leaves in tea plants ( Figure 6) by a mechanism distinct from that previously documented in Arabidopsis (Boachon et al., 2014).
Based on these results, we propose a putative working model for the function of CsUGT87E7 in pathogen infection (Figure 8). When the tea plant is invaded by a fungus, CsUGT87E7 expression is upregulated, leading to the accumulation of SGE. The demand for free SA leads to the upregulation of SA biosynthesis genes in the early disease stages. Consequently, the elevated SA causes rapid ROS accumulation (Figure 7) by inhibiting the activities of ROSscavenging enzymes (Chen et al., 1993). The increased ROS may act as secondary messengers to induce pathogenrelated gene expression (Figure 6), and thus activates the plant disease resistance response (Figure 8). On the other hand, when the expression of CsUGT87E7 is silenced, both SGE and SA formation decreases, dampening the SAassociated resistance response. Taken together, our findings demonstrate that CsUGT87E7 is a SA carboxyl glucosyltransferase that plays a positive role in disease resistance of tea plants by modulating SA homeostasis, although the detailed mechanisms of SGE formation during plant defense require further investigation.

Plant materials and growth conditions
Tea plant samples (C. sinensis var. sinensis cv. Suchazao) were collected from the horticultural research station of Anhui Agricultural University during the early spring. The samples used for metabolite and transcriptome analysis were immediately frozen in liquid nitrogen and stored at -80 C until use.

Gene suppression of CsUGT87E7 in tea plants
Candidate antisense oligonucleotides (AsODN) were selected using SOLIGO software (Ding and Lawrence, 2003) with CsUGT87E7 as the input sequence (Supplemental Table S2). AsODNs were synthesized by the General Biosystems Company (Anhui, China). To silence CsUGT87E7 in leaves Data are presented as the mean ± SD of at least three biological replicates (each with three technical replicates). Significant differences between the treatment and the control group were calculated by one-way ANOVA using SPSS 21.0. *P 5 0.05, **P 5 0.01. that were still attached to the plant, 1 mL of 10 lM AsODN CsUGT87E7 solution was injected into tea seedlings, and seedlings injected with the sODN were used as controls. After 24 h of incubation, the leaves were inoculated as described below.

Pathogen experiments
Identified fungal isolates were cultured on PDA plates at 28 C for 4-5 d. Spores from the culture were suspended in sterile distilled water, the spores were disrupted with a sterile toothpick (Harteveld et al., 2014), and the spore suspension was adjusted to a concentration of 1Â 10E7 conidia per mL using a hemocytometer. Before inoculation, the upper and lower surfaces of the leaves were disinfected with 75% alcohol, rinsed with sterile water, and air-dried. For inoculation, 50 lL of conidial suspension was inoculated onto the upper surface of the tea leaves over three needlepierced holes. Control plants were inoculated with distilled water. Inoculated detached branches were placed in pure water in Erlenmeyer flasks and covered with plastic bags to provide continuous high humidity in the greenhouse (25 C ± 2 C) . To study the systemic disease resistance against Pseudopestalotiopsis, the lower leaves of one-year-old tea plants were primarily inoculated with fungal (Ps. camelliae-sinensis). At 24 h after primary infection, the upper leaves were collected for RNA extraction and SA analysis. Photographs of the disease resistance phenotype of both infected and systemic leaves were taken 4 d after infection of the leaves with pathogens. At least 20 leaves from different plants were assessed for each genotype during each experiment.

RNA isolation, cDNA cloning, and sequence analysis
Total RNA was isolated from leaves of the Shuchazao cultivar using RNAiso-mate for Plant Tissue (Takara, Dalian, China) and RNAiso Plus (Takara) according to the manufacturer's instructions. cDNA was synthesized from total RNA by reverse transcription using HiScript III RT SuperMix for RT-qPCR (Vazyme Biotech Co. Ltd, Nanjing, China). Openreading frame sequences were amplified using Phusion High-Fidelity DNA Polymerase (New England Biolabs, Ipswich, MA). The PCR products were purified with a Gel Extraction Kit (CWBIO, Jiangsu, China), ligated into the pGEX-4T1 vector, and subsequently transformed into Trans1T1 competent cells.

RT-qPCR analysis
RT-qPCR was performed according to published protocols (Jing et al., 2019;Song et al., 2015) using gene-specific primers (Supplemental Table S2). The glyceraldehyde-3-phosphate dehydrogenase (GAPDH) gene was used as an internal reference gene and relative expression was calculated using the 2 -᭝᭝Ct method (Livak and Schmittgen, 2001). All reactions were performed on the CFX96TM System (Bio-Rad, Hercules, CA) using the temperature program 95 C for 3 min and 40 cycles of 95 C for 10 s, 62 C for 30 s.

Heterologous protein expression and purification
The full-length sequence of CsUGT87E7 was digested with BamH1 and Smal1, and the resulting gene fragments were subcloned into the expression vector pGEX-4T-1. The recombinant plasmids were transformed into E. coli BL21 Figure 8 Working model of CsUGT87E7 function in the regulation of SA-mediated defense in tea plant. When the tea plant is invaded by a fungus, CsUGT87E7 expression is up-regulated, leading to the accumulation of SGE. The demand for free SA leads to the up-regulation of SA biosynthesis genes in the early disease stages. As a consequence, levels of pathogen-related proteins increase and plant defense are activated.
(DE3) pLysS cells, which were then cultured at 37 C in Luria-Bertani (LB) liquid medium containing ampicillin (100 lg/mL) and chloramphenicol (25 lg/mL). Isopropyl-b-Dthio-galactopyranoside (IPTG, 1 mM) was added when the optical density (OD600) of the cultured cells reached 0.6-0.8. They were then incubated at 16 C and 150 rpm to permit protein expression. The fusion proteins were purified with GSTÁBind resin (Jing et al., 2019) and protein concentration was determined by a photometric method (Bradford, 1976). The protein sizes were confirmed by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE).

Enzymatic activity assays
In the initial screening, each reaction mixture (25 lL in total) contained Tris-HCl buffer (50 mM, pH 7.5, 10% glycerol [v/v], and 10 mM 2-mercaptoethanol), 250 mM UDPglucose, substrates (20 mM substrate solution), and purified protein (0.5 lg per reaction) as described in Jing et al. (2019). The reaction mixture was incubated at 30 C for 30 min and the reaction was stopped by mixing the reaction solution with the same volume of UDP-Glo assay reagent (Sheikh et al., 2017). The optimum reaction temperature and pH were determined as described in Jing et al. (2019). The optimized conditions were used for subsequent determination of kinetic parameters and at least seven different substrate concentrations and at least three biological replicates were used.

Quantification of SA and SA glucosides
Samples were ground in liquid nitrogen and kept at -80 C before analysis. For metabolite analysis, 50 mg of sample was extracted twice with 1 mL of 75% (v/v) methanol and chlorophenylalanine solution (3 lg/mL) was added as an internal standard. The solution was mixed by vortexing and sonicated at 4 C for 20 min. After that, the mixture was centrifuged at 12,000 rpm and 4 C for 10 min. Finally, the supernatants were collected and used for LC-MS analysis (Jing et al., 2019). Products were identified by comparison of their retention time and MS spectra with those from literature or reference material. The levels of SA, SAG, and SGE were quantified by comparing the signal in the total ion chromatogram with that of the established calibration curve. SA: y = 709.73x + 528.3, R 2 = 0.9988; SAG: y = 149.2x-11036, R 2 = 0.9912; SGE: y = 126.6x -30.85, R 2 = 0.9931.

Determination of ROS content in tea leaves
The hydrogen peroxide (H 2 O 2 ) and superoxide anion (O2 Á-) were determined using hydrogen peroxide (H 2 O 2 ) content assay kit and superoxide anion free radical (O2 Á-) content assay kit purchased from Nanjing Jiancheng Bioengineering Institute (Nanjing, China) according to the manufacturer's instructions.

Statistical analysis
All experiments were carried out with at least three independent biological replicates. Each measurement was carried out in triplicate. Data represent the mean ± SD of three biological replicates. Data were statistically analyzed by one-way analysis of variance (ANOVA) performed using SPSS.

Phylogenetic tree analysis
The phylogenetic tree was constructed by the neighborjoining method using MAGA6 software based on the deduced amino acid sequences of UGTs. The 0.1 scale bar indicates 0.1 nucleotide substitutions per site.

Supplemental data
The following materials are available in the online version of this article.
Supplemental Figure S1. Agarose gel and SDS-PAGE of CsUGT87E7.
Supplemental Figure S2. Glucosyltransferase activity of CsUGT87E7 toward SA with different sugar donors.
Supplemental Figure S3. The effect of different incubation times on the product formation of CsUGT87E7 using SA and UDP-glucose as substrates.
Supplemental Figure S4. The temperature optimization of CsUGT87E7.
Supplemental Figure S5. The pH optimization of CsUGT87E7.
Supplemental Figure S6. The leaf surfaces of tea plants 1 day after fungus infection.
Supplemental Table S1. Screening of enzyme activity toward SA.
Supplemental Table S2. List of primers used in this work.
Conflict of interest statement. None declared.